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IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 

In re Patent Application of 

ENGLAND, J. et al. Atty. Ref.: 1652-13 

Serial No. unknown Group: 
Filed: February 12, 2002 Examiner: 
For: DETERMINING PROPERTIES OF MAGNETIC ELEMENTS 



* * 



February 12, 2002 



Assistant Commissioner for Patents 
Washington, DC 20231 

Sir: 



0 

fy PRELIMINARY AMENDMENT 

M 

jy In order to place the above-identified application in better condition for 

examination, please amend the application as follows: 
IN THE SPECIFICATION 

Please substitute the following paragraphs in the specification for corresponding 
paragraphs previously presented. A copy of the amended specification paragraphs 
showing current revisions is attached. 

Page 1, before the first line, insert as a separate paragraph: 
This application is the US national phase of international application 
PCT/GB00/03092 filed 1 1 August 2000, which designated the US. 
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IN THE CLAIMS 

Please substitute the following amended claims for corresponding claims 
previously presented. A copy of the amended claims showing current revisions is 
attached. 

3. A method according to claim 1, including grouping the transition data by 
the type of element transition. 

6. A method according to claim 4, wherein a signal defining a transition is 
received by one or more receiver coils, including determining the type of transition in 
accordance with the polarity of the rate of change of the field vector in the direction of 
the element. 

7. A method according to claim 4, including determining information relating 
to the switching fields for each of the first and second types of transition. 

10. A method according to claim 8, further comprising calculating the bias field 
on the element as substantially the sum of first and second switching fields. 

11. A method according to claim 9, wherein the first switching field comprises 
the mean value of the switching fields for the first type of transition and the second 
switching field comprises the mean value of the switching fields for the second type of 
transition. 
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12. A method according to claim 1, including associating the transition data 
with one or more respective elements using a receiver vector whose components 
represent the amplitudes of the signals in one or more receive coils. 

13. A method according to claim 1, comprising scanning the tag using a 
rotating magnetic field. 

15. A method according to claim 1, comprising determining the coercivity, the 
local magnetic field bias resolved in the direction of the or each magnetic element and the 
orientation of the or each magnetic element relative to a known interrogation field 
reference frame. 

16. A method according to claim 1, further comprising determining the 
amplitude response of the or each magnetic element to the applied magnetic field. 

19. A method according to claim 17 wherein the respective characteristics 
comprise the coercivities of the elements. 

20. A method according to claim 17, comprising storing data by reference to 
the respective characteristics of the elements. 
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22. A method according to claim 20, wherein data is storable by reference to 
parameters relating to any one or more of rate of change of applied field, perpendicular 
field, response time, characteristic response shape and the statistical distribution of the 
parameters. 

24. A computer program, which when run on a computer, is configured to carry 
out the steps of claim 1. 



the element direction which minimises the scatter of transition point field vectors 



28. A tag reader according to claim 25, wherein the transition data includes 
data defining first and second switching fields at which at an element undergoes first and 
second transitions. 




27. A tag reader according to claim 25, further comprising means for selecting 



0 
U1 

ry 
e 



resolved along the direction of the element. 
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REMARKS 



Attached hereto is a marked-up version of the changes made to the specification 
and claims by the current amendment. The attached page(s) is captioned " Version With 
Markings To Show Changes Made ." 



LSN:ecb 

0 1 100 North Glebe Road, 8th Floor 
Jfl Arlington, VA 22201-4714 
(*f Telephone: (703) 816-4000 
m Facsimile: (703)816-4100 

Hi 



Respectfully submitted, 



a 



NIXON & VANDERHYE P.C. 
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VERSION WITH MARKINGS TO SHOW CHANGES MADE 

IN THE SPECIFICATION 

Page 1, before the first line, insert as a separate paragraph: 

This application is the US national phase of international application 
PCT/GB00/03092 filed 1 1 August 2000, which designated the US. 

U IN THE CLAIMS 

g 

$ 3. A method according to claim 1 or 2 , including grouping the transition data 

I 



m 



8 

CI 



by the type of element transition. 



6. A method according to claim 4-ef-5, wherein a signal defining a transition 



U is received by one or more receiver coils, including determining the type of transition in 
W accordance with the polarity of the rate of change of the field vector in the direction of 
the element. 



7. A method according to claim 4> 5 or 6, including determining information 
relating to the switching fields for each of the first and second types of transition. 

10. A method according to claim 8-of-9, further comprising calculating the bias 
field on the element as substantially the sum of first and second switching fields. 
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11. A method according to claim 9 or 10, wherein the first switching field 
comprises the mean value of the switching fields for the first type of transition and the 
second switching field comprises the mean value of the switching fields for the second 
type of transition. 



ry 

M 

m 



12. A method according to any one of the precedin g claimsj,, including 
associating the transition data with one or more respective elements using a receiver 
vector whose components represent the amplitudes of the signals in one or more receive 
coils. 

13. A method according to any one of "the preceding claims!, comprising 
scanning the tag using a rotating magnetic field. 

15. A method according to any one of the preceding claims J,, comprising 
determining the coercivity, the local magnetic field bias resolved in the direction of the or 
each magnetic element and the orientation of the or each magnetic element relative to a 
known interrogation field reference frame. 



16. A method according to any one of the precedin g claims J[, further 
comprising determining the amplitude response of the or each magnetic element to the 
applied magnetic field. 
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19. A method according to claim 17 or 18 wherein the respective characteristics 
comprise the coercivities of the elements. 

20. A method according to a&y-one of claims- 174e-49, comprising storing data 
by reference to the respective characteristics of the elements. 

22. A method according to claim 20 or 2 1 , wherein data is storable by reference 
to parameters relating to any one or more of rate of change of applied field, perpendicular 
field, response time, characteristic response shape and the statistical distribution of the 
parameters. 

24. A computer program, which when run on a computer, is configured to carry 
out the steps of an y o n e of claims 1 to 23. 

27. A tag reader according to claim 25 or 26 , further comprising means for 
selecting the element direction which minimises the scatter of transition point field 
vectors resolved along the direction of the element. 

28. A tag reader according to aay one of claims 25 to 2 7 , wherein the transition 
data includes data defining first and second switching fields at which at an element 
undergoes first and second transitions. 
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Determining Properties of Magnetic Elements 
Field of the Invention 

This invention relates to magnetic elements, particularly but not exclusively to 
5 methods of distinguishing between magnetic elements and methods and apparatus 
for reading magnetic data tags which include one or more magnetic elements, each 
of which can differ in coercivity, saturated dipole moment (i.e. response amplitude), 
^ orientation and bias field. 

m 

tSfejj 

y 

4* 10 Background 

]g Co-pending PCT publication number WO99/35610 describes tags and reader 

W systems primarily intended for tags fabricated from magnetic material of low 

'T coercivity, with elements at different orientations, in which data is recorded 

primarily by means of the orientation of the elements with respect to each other. 

\ «J ' 15 The described system assumes that the coercivities of the tag elements are all the 

N ■ same, and are very small compared to the interrogation field. 

0 
% 

Summary of the Invention 

According to the invention, there is provided a method of reading a magnetic tag 
20 having at least one magnetic element, comprising interrogating the tag with a 

scanning magnetic field, determining transition data associated with changes in the 
magnetisation state of the at least one magnetic element, associating the transition 
data with one or more respective elements; and for each element, determining the 
element direction which corresponds to the transition data for that element. 



25 



Preferably, the element direction is determined by selecting the direction that 
minimises the scatter of the transition field vectors resolved along the direction of 
the element. 



30 The transition data for each element can be grouped into two sets, which can be 
referred to as forward and reverse transitions. All those in the forward transition 
group have a positive component of the field vector dH/ dt along the element 
vector, and all those in the reverse group have a negative component of dH/dt 
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along the element vector. Mean field values, resolved along the element vector, can 
be calculated. The coercivity of the element is then calculated as half the difference 
between the forward and reverse mean values, and the bias field along the element 
is calculated as the sum of the forward and reverse mean values. 

According to the invention, there is further provided a method of distinguishing 
between a plurality of magnetic elements, comprising the steps of applying a 
scanning magnetic field to the elements, determining the direction of each of the 
elements, for each of the elements, determining the components of the field in the 
direction of the element at which the element switches magnetisation states; and 
from said components, determining, for each of the elements, respective 
characteristics of the element. 

The invention further provides a method of determining, for a magnetic element, 
any one or more of a plurality of characteristics comprising the coercivity of the 
element, the local magnetic field bias resolved in the direction of the element and 
the orientation of the element, comprising the steps of applying a varying magnetic 
field to the element, determining the direction of the element, determining the 
components of the field in the direction of the element at which the element 
switches magnetisation states; and from said components, determining the one or 
more characteristics of the element. 

According to the invention, there is also provided a magnetic tag reader for reading 
a magnetic tag having at least one magnetic element, comprising means for 
interrogating the tag with a scanning magnetic field, means for determining 
transition data associated with changes in the magnetisation state of the at least one 
magnetic element, means for associating the transition data with one or more 
respective elements; and means for determining, for each element, the element 
direction which corresponds to the transition data for that element. 

Brief Description of the Drawings 

Embodiments of the invention will now be described, by way of example, with 
reference to the accompanying drawings, in which: 
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Figure 1 is a schematic diagram of a magnetic data tag reading system; 

Figure 2 is a schematic diagram showing the components of the magnetic data tag 

reading system of Figure 1 in more detail; 

Figure 3 is a schematic diagram showing details of the signal processor/ controller 

illustrated in Figures 1 and 2; 

Figures 4 and 5 illustrate the receive coil set; 

Figures 6 and 7 illustrate the transmit coil set; 

Figure 8 illustrates the antenna comprising transmit and receive coil sets; 
Figure 9 is a flow chart illustrating the overall processing algorithm; 
Figure 10 is a schematic diagram of data acquisition circuitry; 
Figure 1 1 illustrates the transmit current waveforms; 

Figure 12 illustrates a flow chart for the signal processing and filtering algorithm; 
Figure 13 illustrates signals at the inputs to the ADC from the x, y and z receiver 
coil preamplifiers, for a single element transition; 

Figure 14 illustrates the composite filter output of the signal in Figure 13; 

Figure 15 illustrates a 3D scatter plot of the filtered receiver vectors; 

Figure 16 illustrates a flowchart for the clustering algorithm used for planar tags; 

Figure 17 illustrates the composite filter output for three parallel elements; 

Figure 18 illustrates a flowchart for a parallel element clustering algorithm; 

Figure 19 illustrates a 3D scatter plot of the transition field vectors for a single 

element; 

Figure 20 illustrates the same 3D scatter plot as Figure 19, tilted such that the 
transition planes are edge-on; 

Figure 21 illustrates distribution of field vectors that occur along a misaligned 
element direction vector; and 

Figure 22 illustrates a flowchart for the calculation of the mean switching field, 
switching field variance, coercivity and DC bias field. 

Detailed Description 

Referring to Figure 1, a magnetic tag reading system comprises a magnetic data tag 
1, an interrogation unit 2 and a signal processor/controller 3. Magnetic tags 1 to be 
used with a magnetic tag reader according to the invention can record information 
by means of elements of differing coercivities, local bias fields and response 
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amplitudes, as well as orientation. This includes tags described in PCT publication 
number WO99/35610, as well as tags described in, for example, US 5,204,526, US 
5,729,201 and W098/26312. In general terms, magnetic tags 1 comprise magnetic 
elements which typically switch magnetisation state, for example magnetisation 
direction, at given values of applied field depending on element properties, for 
example coercivity. These elements include, for example, thin film elements, 
bistable elements, Barkhausen wire elements and high-permeability elements. The 
applied field which causes switching depends on the magnitude of the component 
of the interrogation field vector in the direction of the element. 

Referring to Figure 2, the tag 1 is attached to an item being labelled or tagged 4, and 
is placed within an interrogation volume 5 within the interrogation unit 2. The 
interrogation unit 2 includes an antenna 6, which comprises transmit and receive 
coil sets 7, 8. The tag 1 is interrogated by a scanning magnetic field 9 generated by 
the transmit coil set 7 under the control 10 of the processor/ controller 3. In 
response to the interrogating magnetic field 9, the tag 1 generates a detectable 
magnetic field response 11, which is detected by the receive coil set 8. The 
processor/controller 3 receives input signals 12, 13 from the transmit and receive 
coil sets 7, 8 respectively and processes the signals to decode data stored on the tag, 
which is made available at an output 14. 

Referring to Figure 3, the processor/ controller 3 comprises a waveform generator 
15 for driving the transmit coil set 7, data acquisition circuitry 16 for receiving 
respective input signals 12, 13 from the transmit and receive coil sets 7, 8 and a 
digital signal processor 17 for processing the resulting output signals 18 from the 
data acquisition circuitry 16 to provide the decoded tag data 14. 

The transmit and receive coil set arrangement 7, 8 is described in detail by reference 
to Figures 4 to 8. 

Figures 4 and 5 illustrate the receive coil set 8. The receiver coils are constructed 
on a cylindrical former 20 of diameter 200 mm and length 400mm. Figure 4 
illustrates the three sets of orthogonal coils used to couple with the tag magnetic 
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elements within the interrogation zone. For the y-direction, the receiver coil set 
comprises 4 coils 21, 22, 23, 24. Inner coils 22, 23 lie on the former 20 and extend 
120mm along the x-direction 25. Both inner coils 22, 23 comprise 100 turns 0.4mm 
ecw. The outer coils 21, 24 comprise 58 turns of 0.4mm ecw and are wound on a 
second co-axial former (not shown) 260mm in diameter. The coils extend 156mm 
along the x-direction. The four coils 21, 22, 23, 24 are connected in series in the 
electrical sense illustrated and 'balanced* by small mechanical re-alignments to 
achieve zero sensitivity to a uniform magnetic field. A second receiver coil set as 
illustrated is sensitive to tag generated field in the z-direction. This coil set is 
identical to the coils 21, 22, 23, 24 but rotated through 90° as shown. The third coil 
set sensitive to tag generated field in the x-direction comprises two solenoid coils 
26, 27. The inner coil 26 comprises 100 turns 0.4mm ecw wound on the former 20, 
and is 120mm long. The outer coil 27 comprises 58 turns of 0.4mm ecw wound on 
the second 260mm diameter co-axial former and is 156mm long. Figure 5 illustrates 
all the coils wound on the inner former 20, and the outer former 28. 

Figure 6 illustrates the three orthogonal transmit coils configuration 7. The coils 
are wound on a cylindrical former 30, 370mm long and 300mm diameter. A 
uniform magnetic field in the y-direction is produced by four coils 31, 32, 33, 34. 
First and third coils 31, 33 comprise a 'modified Helmholz 5 arrangement similar to 
coils 15 and 16. Second and fourth coils 32, 34 comprise a second modified 
'Helmholz' arrangement, with a magnetic axis 25° offset from the first and third 
coils 31, 33. The two 'modified Helmholz' coil sets have magnetic axes 12.5° either 
side of the y-direction. The first coil 31 comprises 50 turns 1.4mm ecw and extends 
370mm in length along the former. Where this coil 31 connects across the open 
end of the former 30, the coil is a flattened half circle with the total coil aperture 
width of 570mm. The two edges of the coil 31 that lie along the solenoid (x- 
direction) subtend 120° at the axial centre of the former. Second to fourth coils 32, 
33, 34 are identical in size and form. Their orientation around the former 30 is 
described above. The four coils are connected in series in the sense illustrated. A 
second transmit coil set generates uniform field in the z-direction. This set 
comprises four identical coils orientated in an orthogonal direction as illustrated. 
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The final transmitter coil consists of a long solenoid coil 35 comprising 260 turns of 
1.4mm ecw on the coil former. This generates uniform field in the x-direction. 

Figure 7 shows the overall transmit coil arrangement for generating uniform field in 
three orthogonal directions. Figure 8 illustrates the antenna 6. The transmit coils 
on the former 30 are located co-axially with the receiver coil tube 20. The 
interrogation volume 5 is defined by a further 190mm ID co-axial tube (not shown) 
that is used to define a mechanical constraint on possible tag positioning in the 
antenna 6. The longitudinal region of highest sensitivity is less than 10cm long and 
tags can be accurately read when separated by 10cm or more along the axis of the 
reader tube. 

Figure 9 illustrates the overall sequence of steps required to decode data stored on a 
magnetic tag. The first stage is data acquisition (step si). Data is acquired by 
detecting the field 11 resulting from the application of a scanning interrogation field 
9 to the tag 1, digitising the resulting signals and storing them for subsequent 
processing. This results in 3 channels of input data, one for each of the x, y and z 
directions. Digital signal processing is carried out to identify individual switching 
points, also referred to herein as transitions (step s2). This results in an array of 
transition information. Each transition is associated with an element (step s3) to 
provide an array of elements. The elements are then individually decoded (step s4). 
Finally, the tag is decoded to provide tag value data (step s5). 

Data Acquisition 

Referring to Figure 10, an example of the signal processor/controller 3 according to 
the invention comprises a National Instruments PCI6711 4-channel DAC card 39 
for waveform generation and a National Instruments PCI6110E 4-channel ADC 
card 40 for data acquisition. The cards are mounted into an industry standard IBM 
compatible PC 41 running Windows 95™. The waveform generation card 39, under 
software control, generates three transmit excitation voltages 42, 43, 44 which are 
passed through respective low-pass filters 45, 46, 47 and amplified by respective 
power amplifiers 48, 49, 50 to drive respective orthogonal transmit coils, which are 
arranged in a series resonant configuration with respective capacitors 51, 52, 53 and 
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resistors 54, 55, 56. The drive current is for example 3A rms to generate a 2.5kA/m 
interrogation field. The transmitter currents are monitored by respective current 
sense resistors 57, 58, 59 which are fed through respective amplifiers 60, 61, 62 as 
inputs 63, 64, 65 to the data acquisition card 40, where they are digitised at a sample 
rate of, for example, 160kHz. The instantaneous transmit field vector can be 
determined from these three signals with knowledge of the relationship between the 
transmit coil field and the current response. For example, pre-calibration of the 
system is carried out by measuring the transmit field for different values of driving 
current. 

Signals induced in the orthogonal receive coils are amplified by respective amplifiers 
66, 67, 68, filtered by respective 130Hz notch filters 69, 70, 71 to remove any 
transmit field component, and fed as inputs 72, 73, 74 to the data acquisition card, 
where they are digitised at a sample rate of, for example, 160kHz. The data 
acquisition is buffered in such a fashion that data is clocked into a buffer, and read 
from the buffer asynchronously at some later point. The buffer depth is sufficient 
to accommodate the worse-case latency in the subsequent processing step. 

A continuous scan is used in this example to interrogate the tag, based on a nominal 
130 Hz rotating magnetic field, whose normal vector is arranged to trace out a spiral 
scan over the surface of a complete sphere, tracing a path from one pole of the 
sphere to the other and back. The equations for the components of the 
'transmitted' B interrogation field are given by: 

B x - (cos 2 (<j>) *cos(9) + sin 2 (<)>)) *cos(e>*) 

+ (sin(<f>)*cos(<j>)*cos(9) - cosfo)*sinfo))*sin(cDt) 

B y - (cos(<j>)*sin(<}>)*cos(9) - sin(+)*cos(*))*cos(cttt) 
+ (sin 2 (<f>)*cos(9) + cos 2 (<(>)) *sin(a>r) 

B z - (-cos(<j>))*sin(9)*cos(cor) - sin(9)*sin(<|>)*sin(a>r) 
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where : is the time, o is the angular frequency of the 130Hz scan, <j> — (constant)* 0, 
and 0 • cos _1 (l-t/7). T is the total time for one complete interrogation. 

Figure 11 shows the three transmit current waveforms 80, 81, 82 received at the 
5 inputs 63, 64, 65 of the data acquisition card 40. 

Digital Signal Processing 

The digital signal processing stage (step s2) performed on the data input to the data 
acquisition card 40 is now described with reference to the flow chart description of 
10 the processing algorithms in Figure 12. 

The purpose of the DSP algorithm is to identify individual transitions, and to record 
all the relevant parameters on each transition for subsequent processing algorithms. 
This leads to a large reduction in the volume of data passed on to the subsequent 
15 processing stages. 

The DSP algorithm operates on the three channels of sample data produced by the 
data acquisition process. Figure 13 shows the raw impulse responses in the x, y and 
z channels for a single transition. 

20 

Referring to Figure 12, in a first step slO, an FIR filter is applied to all three 
channels to produce three sets of filtered data, which form a receiver vector. The 
simplest filter consists of three rectangular sections, and provides a method of 
measuring the height of the peaks in the raw data. If the central section has width 
25 w and height + 1, then the outer two sections have width w/2 and height -1. The 
width, w, is typically the same value as the response time of the magnetic element, 
for example, 20-30us. 

The transmit field vector, H, is used to determine the correct polarity of the element 
30 transition in each of the x, y and z receiver coils. The transition polarity in any 
receiver coil is directly related to the polarity of the rate of change of the field 
vector dH/dt in the direction of the receiver coil. dH/dt values are used to 
produce a "polarity vector", where each component can take the value ±1. The 
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scalar (dot) product of the polarity vector with the filtered receiver vector is 
calculated (step sll) and the receiver vector magnitude, a positive number, is 
multiplied by the sign (±1) of the result (step sl2). This results in the composite 
signal shown in Figure 14, in which the polarity of transitions for every element is 
always the same, allowing the use of a simple peak detector to determine peak 
values. 

Peak detection techniques are well known in the art. In this case, a simple threshold 
is used to gate the peak detector input data, to avoid noise appearing as spurious 
peaks. A peak is identified when three or more values exceed the threshold (step 
sl3), and where the current value is greater than both the previous and next value 
(step si 4). The time of the peak is interpolated to a greater resolution than the 
sample frequency by a simple quadratic fit to these three points. 

The data for each transition is stored in an array (step si 5). The data includes: 
Time 

Field vector (H) 

Rate of change of field vector (dH/dt) 
Receiver vectors (both raw and FIR filtered) 

Element Association 

The function of the element association algorithm is to associate transition data 
points with particular magnetic elements in the tag. Subsequent processing steps 
can then analyse the data for each magnetic element in isolation, thereby reducing 
an apparently complex problem with multiple elements into a series of relatively 
simple numerical solves. 

There are two primary mechanisms that are used to associate transitions with 
elements, depending on whether the elements are generally parallel or not. These 
are described below. In the general case, the first step is to separate into groups 
using a non-parallel algorithm, and then, if required, to analyse each separate group 
to see if it contains more than one parallel or near-parallel elements. 
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For non-parallel elements, the filtered receiver vectors are used to separate out the 
transitions between elements. This can be clearly seen from Figure 15, which shows 
the filtered receiver vectors in a 3D scatter plot for an example tag having seven 
non-parallel elements. Inspection of the plot shows that the majority of the 
transition points lie along one of 7 different lines through the origin, which 
indicates that there are 7 discernible directions of elements in the example tag. 
Each direction can be described by two parameters, and therefore the transitions 
can be clustered together into groups in 2D. There are a number of different 
appropriate techniques than can be used to achieve this multi-dimensional clustering 
(e.g. S. Makeig, S. Enghoff, T-P. Jung, M. Westerfield, J. Townsend, E. Courchesne 
and T.J. Sejnowski, ''Moving- Window Independent Component Analysis of Event- 
Related EEG Data: Component Stability, Journal of Neurophysiology"). Additional 
knowledge about the particular tag construction can be useful to simplify the 
problem. For example, if all the elements are in the same plane, then the problem 
can be reduced to a one-dimensional problem. Knowledge of the number of 
elements expected can assist in the clustering process. 

In the particular case of a planar tag, with a known number of elements, the 
algorithm outlined in Figure 16 is used. The normal to the plane of the transitions 
20 is determined by, for example, a numerical process (step s20). For example, the dot 
product of every receiver vector with an estimated direction vector is calculated, 
and this process is iterated until the sum of the magnitude of the dot products is 
minimised. This reduces the problem to a 1-D problem i.e. the angle in the plane. 
An in-plane set of vectors can be calculated from the original set of vectors simply 
25 by subtracting from each vector in turn the dot product of itself with the normal to 
the plane. In-plane angles between any two in-plane vectors can them simply be 
calculated in the usual way using dot products. All the in-plane angles are wrapped 
into the range 0-180° by adding or subtracting multiples of 180° as required. The 
algorithm calculates a histogram of in-plane angles relative to some arbitrary datum, 
30 such as the first point. For example, if the histogram bins are 1° wide, then the nth 
bin will contain a count of the number of angles that fall in the range n° to (n+ 1)°. 
This will typically give a series of peaks, one for each element. For example, the 
algorithm obtains the second point from the transition array (step $21), measures 
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the in-plane angle relative to the first point (step s22) and increments the 
appropriate bin of the histogram (step s23). This process is repeated until all the 
data has been processed (step s24). After applying Gaussian smoothing to the 
histogram data (step s25), the direction of an element in the tag can be found by 
determining the highest peak in this histogram (step s26). To determine the 
transitions that belong to the element in a given direction, the algorithm finds all 
the transitions that are within, say, 2° (in plane) of this direction (step s27). The 
processed peak is then eliminated from the calculation (step s28) and this 
processing sequence is repeated (steps s26 to s28) until all the data has been 
processed (step s29). 

To separate parallel elements, the algorithm makes use of two properties of a 
continuous scan of the field vector, H, around the elements, first that the elements 
transition in order from the lowest to the highest coercivity and second, that the 
field vector, H, rotates by at least 90° between the last transition of the highest 
coercivity element in one direction and the first transition of the lowest coervicity 
element in the reverse direction. 

If some of the elements do not change state, because the transmit field does not 
reach a high enough value, then there will be fewer transitions, in a 180° scan, than 
there are elements. In this case, transitions are "lost", starting with the highest- 
coercivity element. Figure 17 shows the filtered composite waveform over a few 
rotations of transmit field, for three parallel elements with different coercivities. 
Each element is associated with a respective peak 90, 91, 92 and it is relatively 
straightforward to separate out the transitions belonging to different elements. An 
outline algorithm to achieve this is shown in Figure 18. 

This works by maintaining an element counter that is incremented each time a new 
transition is identified, and set to zero each time the field rotates by more than 90° 
between transitions. Data is extracted from the transition array (step s30) and the 
algorithm determines whether the transmit field has rotated by more than 90° since 
the last transition point (step s31). If it has, the element counter is reset to zero 
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(step s32). Following this, the element number corresponding to the transition is 
set according to the current value of the element counter (step s33). The element 
counter is then incremented (step s34) and the process repeated for the next 
transition point (step s30). For example, the first transition following the zeroing of 
the element counter is associated with element 0, the next with element 1 and so on, 
until the field has rotated by more than 90 degrees. 

Element Decode 

The purpose of the element decoding algorithm is to take transition data belonging 
to one element, and to determine the best-fit direction vector for this element. 
Once the direction is known, the coercivity of the element, and any net DC field or 
"bias" along the element vector can be calculated. 

Figure 19 is a 3D scatter plot of transition point field vectors, for a single bistable 
magnetic element with a finite coercivity. In this example, the field has been 
scanned approximately over the surface of a sphere, so the transition points lie 
roughly on two circles 93, 94. More generally, the transition points would be 
expected to lie on one of two planes. By tilting the view of the transition data in 
the scatter plot, it is possible to show the two planes edge-on, as illustrated in 
Figure 20. The bold vertical arrow 95 shows the element vector. 

The element decoding algorithm attempts to determine the best vector direction for 
the element, by minimising the scatter of field vectors resolved in this direction. 
Figure 21 illustrates the situation where a guess 96 has been taken for the element 
vector that is not in the correct direction. Taking the upper set of transitions 93, it 
is clear that when these are projected onto the element vector 96, they form an 
extended distribution 97 (shown by a darkened section) along the vector 96. As the 
vector is rotated around, the extent of this distribution will be smallest when the 
vector is closest to the actual direction of the tag element. 

Figure 22 shows a flowchart for the algorithm used to calculate the error used in the 
iterative solving process. The algorithm uses the current guess for the element 
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veaor direction, V. Initially, this is the vector direction from the element 
association algorithm. 

The first data point is retrieved (step s40) and the dot product of the direction 
estimate V with the polarity vector for dH/dt is calculated, as described in relation 
to Figure 12 above (step s41). If this is positive (step s42) i.e. for the upper set of 
transitions, then the dot product of the field veaor and the element veaor is 
calculated (step s43) and added to the upper set of statistics (step s44). The dot 
0 produa resolves the component of the field veaor along the element veaor. The 

% 10 same calculation is carried out for the lower set of transitions (steps s45, s46), 

*$* 

# indicated by the negative dot produa at step s41. The upper and lower (forward 

and reverse) sets of transitions are distinguished by the sign of dH/ dt along the 
direaion of the element, or alternatively by the sign of the filtered receiver veaor 
along the direaion of the element. This procedure is repeated for all the data 
Iff j j points (step s47). An average value of variance is calculated from the variances for 

each of the upper and lower sets of transitions separately, weighted by the number 
Q of transitions in the upper and lower sets of transitions (step s48). A standard 

formula is used to calculate the variance for each set of data. For a set of measured 
data points, x, the variance, var(x) is the mean of the squares of x minus the square 
20 of the mean of x, or mathematically 

var(*) = (**)-(x) 2 

The weighted variance of N u upper transition points, «, with variance var(n) and 
25 N t lower transition points, /, with variance var(/) is then given by: 

N u v ar(u) + N f var(7) 



fy 



'4 
rij 



30 



The weighted variance is used as a measure of the error in the guessed veaor 
direaion. When the guessed direaion is equal to the actual element direaion, the 
weighted variance will generally have its minimum value. The value will never fall 
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to zero, because there is always a certain amount of noise in the determination of 
the transition field, arising from sources such as electronic noise and randomness in 
the material behaviour. In the simplest case, the variance is a function (numerically 
evaluated, rather than an analytic function) of two direction variables, such as 8 and 
$ from the spherical polar co-ordinates (r, 9, The value of this function can be 
minimised using a standard numerical minimisation algorithm. The variance varies 
approximately quadratically with the deviation from the ideal direction, and this 
means that the minimisation algorithm can be extremely efficient (the a quadratic" 
case is generally considered to be the easiest). Muiti-variate numerical minimisation 
algorithms are well known in the art - for example, Powell's method. 

When the weighted variance is not minimum (step s49), the direction estimate V is 
adjusted in accordance with the appropriate minimisation algorithm (step s50) and 
the algorithm re-run with the new value of V. When the weighted variance is 
minimised, the mean values of the field for the upper and lower sets of transitions 
are calculated (steps s51, s52). The coercivity of the element is calculated as half 
the difference between the two switching fields (step s53), while the DC field along 
the element is calculated as the sum of the two switching fields (step s54). 

Additional parameters as well as direction can usefully be added to the numerical 
minimisation. The most important term to add is the vector velocity, which allows 
the algorithm to deal with the movement of the tag elements during the decoding 
process. The element direction in steps s41, s43, s45 is then a function of the time 
at which the transition occurs, and the function for the variance then depends on 
four parameters (for example 9, <|>, d9/dt and d<|>/dt). Once again, this function can 
simply be minimised using a standard multi-variate minimisation algorithm. 

Many magnetic elements do not behave ideally, and show significant changes in 
their switching field (or coercivity) depending on the value of dH/ dt. A general 
form for the switching field, H swltch , resolved along the element is: 

„ , (dHY . (dHY 
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where a,b etc are arbitrary powers. If the coefficients k are known, then the value 
of Ho can be calculated from the measured switching field, and the variance of this 
value can be minimised as before. If the coefficients k are not known, but the 
values a, b are known, then the function for the numerical minimisation can also 
include the coefficients, k, as function arguments, as well as the direction and 
velocity terms. In this case, the values of k can be used to distinguish between 
different types of materials and thereby store more data. 

Anisotropic thin-film magnetic materials can exhibit a further form of non-ideal 
behaviour. For materials with an easy-axis of magnetisation, the in-plane field 
perpendicular to the easy axis can influence the field at which the material switches. 
A similar approach to the one described above can be used to calculate a nominally 
constant value, Hq, from the raw switching points. 

Tag Decode 

The primary output data for each magnetic element from a reader according to the 
invention is as follows: 

Orientation in the reader (vector) 
Coercivity of the element (scalar) 
Bias field along each element (scalar) 
Amplitude response (scalar) 

Secondary data for each element includes 

dH/dt coefficients 
Perpendicular field coefficients 
Response time 

Characteristic response "shape" or spectrum 
Statistical distribution of primary parameters 

This data assumes little about the construction of the tag. The structure of the tag 
(e.g. which elements share bias magnet elements) may be used to provide more 
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detail - for example, the magnitude and direction of an overall bias field. The 
details of the chosen coding scheme are used to translate all these raw parameters 
into useful data stored on the tag. 

The above examples of the invention are intended to be illustrative, rather than 
restrictive. A person skilled in the art would understand that various modifications 
and variations in the detailed implementation are possible, and are considered to be 
within the scope and spirit of the invention as defined in the appended claims. 
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Claims 

1. A method of reading a magnetic tag having at least one magnetic element, 
comprising: 

interrogating the tag with a scanning magnetic field; 

determining transition data associated with changes in the magnetisation 

state of the at least one magnetic element; 

associating the transition data with one or more respective elements; and 
for each element, determining the element direction which corresponds to 

the transition data for that element. 

2. A method according to claim 1, wherein the step of determining the element 
direction comprises selecting the direction which minimises the scatter of transition 
field vectors resolved along the direction of the element. 

3. A method according to claim 1 or 2, including grouping the transition data 
by the type of element transition. 

4. A method according to claim 3, comprising grouping first and second types 
of element transition. 

5. A method according to claim 4, wherein the first type of element transition 
comprises a forward transition and the second type of element transition comprises 
a reverse transition. 

6. A method according to claim 4 or 5, wherein a signal defining a transition is 
received by one or more receiver coils, including determining the type of 
transition in accordance with the polarity of the rate of change of the field vector in 
the direction of the element. 

7. A method according to claim 4, 5 or 6, including determining information 
relating to the switching fields for each of the first and second types of transition. 
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8. A method according to claim 7, comprising determining element 
characteristics from said switching field information relating to transition data 
associated with an element. 

9. A method according to claim 8, further comprising calculating the coercivity 
of the element as substantially half the difference between first and second 
switching fields. 

10. A method according to claim 8 or 9, further comprising calculating the bias 
field on the element as substantially the sum of first and second switching fields. 

11. A method according to claim 9 or 10, wherein the first switching field 
comprises the mean value of the switching fields for the first type of transition and 
the second switching field comprises the mean value of the switching fields for the 
second type of transition. 

12. A method according to any one of the preceding claims, including 
associating the transition data with one or more respective elements using a receiver 
vector whose components represent the amplitudes of the signals in one or more 
receive coils. 

13. A method according to any one of the preceding claims, comprising scanning 
the tag using a rotating magnetic field. 

14. A method according to claim 13, in which the tag comprises a plurality of 
magnetic elements, further comprising associating transition data with respective 
elements in accordance with the order in which the elements transition in response 
to the rotating field. 

15. A method according to any one of the preceding claims, comprising 
determining the coercivity, the local magnetic field bias resolved in the direction of 
the or each magnetic element and the orientation of the or each magnetic element 
relative to a known interrogation field reference frame. 
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16. A method according to any one of the preceding claims, further comprising 
determining the amplitude response of the or each magnetic element to the applied 
magnetic field. 

17. A method of distinguishing between a plurality of magnetic elements, 
comprising the steps of: 

applying a scanning magnetic field to the elements; 
determining the direction of each of the elements; 

for each of the elements, determining the components of the field in the 
direction of the element at which the element switches magnetisation states; and 

from said components, determining, for each of the elements, respective 
characteristics of the element. 

18. A method according to claim 17, comprising determining first and second 
switching components as the components when the rate of change of the field along 
the direction of the element is positive and negative respectively. 

19. A method according to claim 17 or 18 wherein the respective characteristics 
comprise the coercivities of the elements. 

20. A method according to any one of claims 17 to 19, comprising storing data 
by reference to the respective characteristics of the elements. 

21. A method according to claim 20, wherein data is storable by reference to any 
one or more of orientation of the elements, coercivity, bias field along the element 
and amplitude response. 

22. A method according to claim 20 or 21, wherein data is storable by reference 
to parameters relating to any one or more of rate of change of applied field, 
perpendicular field, response time, characteristic response shape and the statistical 
distribution of the parameters. 
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23. A method of determining, for a magnetic element, any one or more of a 
plurality of characteristics comprising the coercivity of the element, the local 
magnetic field bias resolved in the direction of the element and the orientation of 
the element, comprising the steps of: 

5 applying a varying magnetic field to the element; 

determining the direction of the element; 

determining the components of the field in the direction of the element at 
which the element switches magnetisation states; and 

from said components, determining the one or more characteristics of the 
10 element* 

24. A computer program, which when run on a computer, is configured to cany 
out the steps of any one of claims 1 to 23. 

15 25. A magnetic tag reader for reading a magnetic tag having at least one 
magnetic element, comprising: 

means for interrogating the tag with a scanning magnetic field; 
means for determining transition data associated with changes in the 
magnetisation state of the at least one magnetic element; 
20 means for associating the transition data with one or more respective 

elements; and 

means for determining, for each element, the element direction which 
corresponds to the transition data for that element. 

25 26. A tag reader according to claim 25, wherein the scanning field comprises a 
rotating magnetic field. 



30 



27. A tag reader according to claim 25 or 26, further comprising means for 
selecting the element direction which minimises the scatter of transition point field 
vectors resolved along the direction of the element. 
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28. A tag reader according to any one of claims 25 to 27, wherein the transition 
data includes data defining first and second switching fields at which at an element 
undergoes first and second transitions. 
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(57) Abstract: A method of determining differing characteristics 
of magnetic dipole elements such as orientation, coercivity, bias 
and response amplitude and a tag reader for reading magnetic tags 
containing such elements. The elements are scanned by a rotating 
magnetic field and two sets of transition data are determined. The 
transition data sets are associated with respective elements and analysed 
to determine mean field values resolved along the element vectors. 
These field values are used to determine properties of the elements, 
such as coercivity. 
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Nixon & Vanderhye P.C. (10/99) 
(Domestic Non-Assigned/Foreign) 

RULE 63 (37 C.F.R. 1.63) 
DECLARATION AND POWER OF ATTORNEY 
FOR PATENT APPLICATION 
IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 

As a below named inventor, I hereby declare that my residence, post office address and citizenship are as stated below next to my name and I believe 
I am the original, first and sole inventor (if only one name is listed below) or an original, first and joint inventor (if plural names are listed below) of the 
subject matter which is claimed and for which a patent is sought on the invention entitled: 
DETERMINING PROPERTIES OF MAGNETIC ELEMENTS 

the specification of which (check applicable box(s)): ~~ "~ " 

□ is attached hereto 

g was filed on as U.S. Application Serial No. 10/049,420 (AttyDkt. No ) 

□ was filed as PCT International application No. PCT/GBOO/03092 on 11 Au gust 2000 ~~ — 

and (if applicable to U.S. or PCT application) was amended on 



I hereby state that I have reviewed and understand the contents of the above identified specification, including the claims, as amended by any 
amendment referred to above. I acknowledge the duty to disclose information which is material to the patentability of this application in accordance 
Mh 37 C.F.R. 1.56. I hereby claim foreign priority benefits under 35 U.S.C. 1 19/365 of any foreign application(s) for patent or inventor's certificate 
^sjed below and have also identified below any foreign application for patent or inventor's certificate having a filing date before that of the application on 
which pnonty is claimed or, if no priority is claimed, before the filing date of this application- " 
Wlofity Foreign Apptrcatidn(s): _ _______ 

l P S^ mber C <r*ry Day/MonthTYear Filed 

*£ gglgigg-g OB 12 August 1999 

thereby claim the benefit under 35 U.S.C. §11 9(e) of any United States provisional application(s) listed below 
Application Number Date/Month/Year Filed 

— iir-r — — . _ 



thereby claim the benefit under 35 U.S.C. 120/365 of all prior United States and PCT international applications listed above or below and insofar as 
f* Va J , 6 ^ ? ?' f I _ ! C ' aimS ° f this a PP |ication is not disclosed in such prior applications in the manner provided by the first paragraph of 35 
UiS.C. 1 12, 1 acknowledge the duty to disclose material information as defined in 37 C.F.R. 1 .56 which occurred between the filing date of the prior 
applications and the national or PCT international filing date of this application: 

? Up. 

lMorU.S./PCTApplication(s): «...„_. . . . 

plication Serial No. Day/Month/Year Filed pending ab^on^ 



I hereby declare that all statements made herein of my own knowledge are true and that all statements made on information and belief are believed to 
be true; and further that these statements were made with the knowledge that willful false statements and the like so made are punishable by fine or 
imprisonment, or both, under Section 1001 of Title 18 of the United States Code and that such willful false statements may jeopardize the validity of the 

ss ^oir A %™z vi z^rr-, A " d on beha i! of the ° wner < s > here ° f - » app°iflLHiaauAj_M___^^ R d . ( 

8 Floor Arimgton VA 22201-4714, telephone number (703) 816-4000 (to whom all communications are to be directed , and rnifoTlowing— 
attorneys thereof (of the sameaddress) individually and collectively owner's/owners' attorneys to prosecute this application and to transact all business 
(m the Paten and Trademark Office connected therewith and with the resulting patent Larry S. Nixon___5S42_Arthur R. Crawford „2532*James T 
I of S T er ?21S£Bnbert W. Fans, 31352; Rjchard G. Besha, 22770; Mark E. Nusbaum,____34___Mich4el J tKeenan, 32106- Bryan H Davidson 30251- 

I Stanley C Spooner pc^onard C. Mitchard, 2__.ca____.ane M. Byers. ^smJeffryt^son, 3X>4__l__Mn R- LasTova 5__iS WarSrf^ 

Burnam Jr^JMary J. Wlson v 3_^_55__J. Scott Davidson , 33489: Alan M. Kagen___m_^ RobertT^!.^. J. bido ^663 J ame S D 
Berqu«st, 3_____60lpdeep S. Gill, 3Z334^4ichael J. Shea^______&.Donald L Jackson, «_2S_U.ichelle N. Lester__SfFrank P P^_SlS Joseoh 

55££5£££ A „*^^**l™°* * ^ ^ris Com_,tzis. ^7Toa ,v T. Tanigaw^cl I also aSS5^ 
Vanderhye to delete any attorney names/numbers no longer with the firm and to act and rely solely on instructionsTJifecfly communicated from the 
person, assignee, attorney, firm, or other organization sending instructions to Nixon & Vanderhye on behalf of the owner(s). 

1. Inventor's Signature: A/^\ fk^7 „ , A Date: ff /f /0Z 

. lnventor: _ J a^es Mark. Carson ..ENGLAND Britiih 

1^0 , , Camhfg ■■ (state/country, fig. Kingdom Q ^ dto "^> 

Post Office Address: 44A Butt Lane, Milton 2 W ^ 

(Zip Code) CB4 6DG 

2. Inventor's Signature: Date- 

,nventor: Ar, drew Nicholas DAMES ' British 

Residence: (city) Cambridge^ ™ (state/countryS' ^United Kingdom Citizenship) 

Post Office Address: 74 De FreviHe Avenue ■ — 



(Zip Code) CB4 1HU 



FOR ADDITIONAL INVENTORS, check box □ and attach sheet with same information and signature and date for each. 
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with 37 C.F.R. 1.56. I hereby claim foreign pnonty benefits under 35 U.S.C. 1 19/365 of any foreign applications) for patent or inventor's certificate 
Ss^ed below and have also identified below any foreign application for patent or inventor's certificate having a filing date before that of the application on 
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iftority Foreign Application^): 
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Lhereby claim the benefit under 35 U.S.C. 120/365 of all prior United States and PCT international applications listed above or below and insofar as 
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tlication Serial No. Day/Month/Year Filed pJJS^SSS 



L^n^ y fn^,^ LJ statemente made here.n of my own knowledge are true and that all statements made on information and belief are believed to 
be true, and further that these statements were made with the knowledge that willful false statements and the like so made are punishable by fine or 
imprisonment, or both, under Section 1001 of Title 18 of the United States Code and that such willful false statements may jeopanlfce th fvalidih of the 
application or any patent issued thereon. And on behalf of the owner(s) hereof, I hereby appoint NIXON & VANDERHYE P C 1100 I North Glebe Rd 
T"l 0r - ?"' f )f A ™-«JK telephone number (703) 816^000 (to whom all communications arc ,tol£ SS), aXEetS? " 

attorneys thereof (of the same address) individually and collectively owner's/owners' attorneys to prosecute this application and to transact allbusiness 
in the Patent and Trademadc Office connected therewith and with the resulting patent Larry S. Nixon, 25640; Arthu ^R CrTXd 25^T?ames T 
Hosmer 30184; Robert W Faris, 31352; Richard G. Besha, 22770; Mark E. Nusbaum, 32348; Michel J. Keenan, 32106 Biyari H D^viS 30251 ■ 
Stanley C Spooner 27393; Leonard C. Mitchard, 29009; Duane M. Byers. 33363; Jeffry H. Nelson, 30481; John R. Lastova 33149 H Warren 
Burnam Jr 29366; Mary J. Wilson, 32955; J. Scott Davidson, 33489; Alan M. Kagen, 36178; Robert A. Molan, 29834; B J Sadoff 36663 jlmes D 
Berquist, 34776; Updeep S. Gill, 37334; Michael J. Shea, 34725; Donald L. Jackson, 41090; Michelle N Lester 3233 Frank I P Presta 1982^1^ 
S. Presta, 35329; Joseph A. Rhoa. 37515; Raymond Y. Mah, 41426; Chris Comuntzis. 31097; Gary T. I abo auTho^l N fxo'n f P 

Vanderhye to delete any attorney names/numbers no longer with the firm and to act and rely solely on instructions direcUy communicSed from the 
person, assignee, attorney, firm, or other organization sending instructions to Nixon & Vanderhye on behalf of the 3) 

1 . Inventor's Signature: 



. , — . Date: 

lnventor: Jam es Mark. Carson ENGLAND British 

Residence- wt«^ r^mhrirtl^ Ml (,ast > (citizenship) 

Res dence. (aty) Cambndge (state/country) United Kingdom 

Post Office Address: 44A Butt Lane, Milton ~ ~ 



(Zip Code) CB4 6DG, V — 

2 inventor's Signature: Date: tSjf^UTW ? 
jW>lnventor: ■ r _ P AMFS Sritiih ~ 



£AM£2 British 



Yfirsfl — — — ~ llir " - w onusn 

— ce < cit > > Cambridge (XMcoun ^l^ Kingdom QS)^^ 

74 He Frm/illo Auonna - ■ — ■ ■- 



Post Office Address: 74 De Freville Avenue 
(Zip Code) CB4 1HU 



FOR ADDITIONAL INVENTORS, check box □ and attach sheet with same information and signature and date for each. 



